Localize energy in random media: A new phase state 
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This Letter reports a new wave phenomenon. We show that wave energy can be stored in random 
media. Associated with such energy storage is a global collective behavior. The features depicted 
are so general that they may be observed in many random systems. 



PACS numbers: 43.20., 71.55J, 03.40K (Chin. J. Phys, 2000) 



A major task in physics research is to search for com- 
mon principles behind various nature processes. A num- 
ber of greatest advances in the last century include un- 
derstanding of how the manifestation of many body in- 
teractions may lead to several important ubiquitous phe- 
nomena including superconductivity, quantum Hall ef- 
fect, screening effect in plasma, and the Kondo effects [|]. 
It becomes a daily experience that macroscopic objects 
arise from interaction between individuals of microscopic 
scale in a many body system . 

In this Letter, we will show a new general wave phe- 
nomenon in certain many body systems. Namely, wave 
energy can be stored in a number of different disordered 
systems. When this occurs, a new phase state emerges. 

Up to date, the fundamental equation describing an 
arbitrary many body system can be more or less written 
as 
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in which £ is called the divisor g , Fi denotes the external 
stimulation, is an operator that can be of temporal, 
spatial or both nature, depicting the interaction between 
consisting bodies. In the linear approximation, the oper- 
ator T is represented by a wave propagator. For example, 
in the Yang-Feldman representation of field equations, T 
is the usual ferminic or bosonic two-point Green's func- 
tion. In Eq. (^, i is the index of the individuals in the 
system of total N bodies. 

So far, there are two known types of divisor [||. Type I 
divisor usually refers to fermionic systems and is written 
as 
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examples. Eq. (|l]) can also describe multiple scattering of 
classical waves such as acoustic propagation in air-filled 
bubbles Q and cylinders and electromagnetic waves in 
electric dipolar systems M. In these systems, the inter- 
action between constituents is mediated by waves which 
themselves are described by Type II equations. 

Here we show some new generic features with Type 
II systems. Consider wave propagation in a Type II 
system. The propagating wave, which can be either 
acoustic or electromangetic, will be interfered by scat- 
tered waves from each constituent being stimulated by 
the incident wave. The energy flow in the system is 
J ~ Re['0(-i)V-0]. Writing the field as V = iV'le*'', the 
current becomes |'0pV6', a version of Meissner equation. 
It is clear that when 9 is constant while ^ 0, the flow 
stops and the energy must be localized in space, hinting 
at that a phase transition occurs and modes are con- 
densed in the real space. This implies that energy can 
be stored in certain spatial domains. Clearly, the con- 
stant phase 9 indicates the appearance of a long range 
ordering in the system. We note that in the supercon- 
ductivity case, the uniform phase (V0 — 0) and a finite 
super-current lead to an infinite density of carriers, indi- 
cating the condensation of electron pairs in the momen- 
tum space and resulting in the infinite conductivity. 

The above simple argument naturally leads to the 
question of whether there are systems which can real- 
ize the above features. We found that these features can 
appear in a number of different systems. The previously 
considered acoustic propagation in bubbly water |^ is 
one of such systems, as well as acoustic propagation in 
water with air-filled cylinders 1^, and the electric diople 
system. 

For all the systems mentioned, in the linear approxi- 
mation the governing equation can be reduced to 



Type II divisor refers to boson-like systems and is written 
as 
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In these, e(V) and aj(V) are the operators yielding the 
energy spectra. 

Equation (^ accommodates a variety of problems of 
great interest. To name a few, this includes the vibration 
of the lattice formed by atoms, spin arrays, and electrons 
in lattices. Interested readers may refer to for more 
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where Fi represents the external stimulation, C is the 
coupling constant, d denotes the dimension, and c is the 
wave speed. Strictly, for the 2D {d = 2) acoustic case, the 
propagator on the RHS would often be in a form of the 
zero-th order Hankel function of the first kind ||^. The 
quantity refers to either the oscillation of the dipoles 
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in electric systems and or to the vibration of the air- 
filled bodies in the acoustic situation, and luq refers to 
the corresponding natural frequency, and 7 the possible 
damping effects due to such as radiation. The parameters 
will depend on specific models considered. For acoustic 
scattering in bubbly water, for example, C ^ a, 7 ~ wq, 
and LO ^ ^y/Po- Here a is the radius of air-bubbles and 
po is the ambient pressure. In this case, ip refers to the 
radial pulsating of the bubbles. The parameters can be 
found in [Q and references therein. For the electric dipole 
systems, C = "f?™ '^i^h q being the charge and m the 
mass of the charge, and V' is the electric dipole. 

Eq. reveals that when it is excited by the incoming 
wave, either of electronmagnetic or of acoustic nature, 
the constituent will radiate or scatter waves. The radi- 
ated or scattered waves will be again re-radiated or scat- 
tered by other constituents, a process repeated to estab- 
lish an infinite recursive pattern of multiple interactions. 
We note here that Eq. (|) can in fact also describe many 
other systems not being considered here. 

To explore the general properties of the model in 
Eq. (^). We solve it numerically in the frequency do- 
main. The generic set up is as follows. There are N 
bodies randomly distributed in a system, and they are 
excited by a point source located roughly in the middle 
of these bodies. We further assume that all elements in 
the system are identical. We write i/^i = Aie^^\ where we 
dropped the time factor e~*"*. 

For each phase 9i , by analogy with [|j we define a phase 
vector such that Vi — cosOiCx + sin^^ej,. The phase vec- 
tors would indicate the degree the coherence of the vi- 
bration behavior. We have done a few general numerical 
calculations for Eq. (^. The significant discoveries we 
found can be summarized as follows. For two or three di- 
mensions, the energy is localized near the source within a 
range of frequencies slightly above the natural frequency 
for sufficiently large coupling constants and number den- 
sities of the constituents. When the energy is localized in 
the system, all the constituents oscillates in phase, i. e. 
the phase 9 becomes constant, revealing a long range or- 
dering in the system. 

Here we would like to show one example. Randomly 
put N identical electric dipoles on the x — y plane, inside 
a box of side length L. All the dipoles points to the pos- 
itive z-axis. Assume that the averaged distance between 
dipoles is d; thus the number density of the dipoles is 
n ~ l/dP. In the computation, the damping factor 7 can 
be adjusted to reflect various situations. A transmitting 
source is put at the center of the dipole cloud. As this 
is the first step in our research, we simplify further that 
all the dipoles can only oscillate along their axes. The 
dipoles will oscillate vertically in response to the source 
and as well as the radiated waves from other dipoles, i. e. 
all the dipoles are coupled via the radiated waves. In 
the computation, all frequencies are scaled by the nat- 
ural frequency of the dipole. The coupling constant C 
and damping 7 can be varied. We study the behavior of 
the aforementioned phase vectors and the energy distri- 
bution in the system; it is easy to see that the energy 



can be represented by the oscillation amplitude of the 
dipoles. 

Fig. |l| illustrates our findings for one arbitrary random 
configuration. For low frequencies, e. g. at uj/liOo = 0.98, 
the phase vectors points randomly. The distribution of 
the energy spreads. When increasing the frequency to 
certain values slightly above the natural frequency, all 
phase vectors tend to point to a uniform direction, in- 
dicating a new phase state. Meanwhile the energy is 
centered at the site of the transmitting source. This is 
shown by the case at oj/ujq = 1.1. When the frequency is 
increased further, the phase ordering disappears, and the 
energy distribution becomes extended again, as shown by 
the case with lo/loq = 10. These features are valid for a 
range of the coupling constant C and the damping rate. 
The parameters used in the present specific example are: 
C = 0.1, n = 0.032, and 7 = 0.0001. We take N as 
900; in Fig. 1, however, only a fraction of the system is 
shown in order to expose the features in the most explicit 
way. The features tend to diminish as either the couple 
constant or the number density is reduced to a certain 
value. When the damping rate 7 is increased to a large 
value, the phase ordering will then be degraded. The 
phase ordering in the present system has a clear physi- 
cal meaning. That is, it refiects that all elements in the 
many body system oscillate or vibrate coherently. We 
stress that for simplicity we have choosen the 2D arrange- 
ment of the dipoles in a 3D system. Strictly speaking, 
due to radiation into the third dimension, waves cannot 
be stored permanently in the medium. The energy will 
gradually escape. We have also performed calculation for 
a 3D random configuration of the dipoles. The similar 
results hold. 

In summary, we have exhibited a new phase state pos- 
sibly observable in a class of many body systems. It is 
shown that energy can be stored coherently in random 
media. The discovery reported here may provide insight 
to the long standing problem of the Anderson localization 
of waves in disordered media. 
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